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1.  Introduction 


This  paper  pertains  to  a  device  for  generating  useful  elec¬ 
trical  power  by  means  of  an  electrohydrodynamic  (EHD)  process. 

The  device  utilizes  the  flow  of  a  fluid,  normally  a  gas,  such  as 
air  or  water  vapor,  in  which  are  entrained  a  very  large  number  of 
very  fine  and  well  distributed  solid  or  liquid  particles,  for 
example,  water  droplets  as  in  an  aerosol  spray.  The  particles  are 
electrically  charged,  either  positively  or  negatively,  as  may  be 
convenient  in  a  particular  case.  The  gas  is  caused  to  flow  through 
a  nozzle-like  channel  by  the  imposition  of  a  suitable  pressure  drop. 
Let  station  0  be  the  inlet  to  the  nozzle.  The  charged  particles 
are  introduced  into  the  stream  of  gas  by  a  suitable  injector  at 
some  upstream  location,  call  it  station  i,  and  are  removed  from  the 
stream  by  a  suitable  collector  grid  at  some  downstream  location, 
call  it  station  2.  An  important  intermediate  location,  station  1, 
is  at  the  throat  of  the  nozzle.  The  electrical  charges  which  are 
transported  downstream  along  with  the  particles  on  which  they  re¬ 
side  constitute  an  electrical  current.  These  charges  move  through 
an  electrical  field  which  exerts  forces  upon  them  in  a  direction 
and  sense  opposed  to  the  general  fluid  motion.  By  proper  design  of 
the  injector,  the  size  of  the  water  droplets  or  particles  may  be  so 
regulated  that  these  droplets  have  low  mobility,  that  is,  they  move 
approximately  with  the  surrounding  gas  with  negligible  relative  slip 
produced  by  the  opposing  electrical  forces. 

The  gas  stream  does  work  on  the  charged  particles  in  moving 
them  against  the  resistance  of  the  electrical  forces.  In  this 
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process  the  gas  stream  undergoes  a  corresponding  decrease  in 
total  enthalpy.  The  work  done  upon  the  charged  particles  creates 
a  difference  in  electrical  potential  between  stations  i  and  2. 

These  stations  are  connected  by  an  external  circuit  which  includes 
a  useful  electrical  load.  Thus  the  enthalpy  drop  of  the  gas  is 
ultimately  converted  into  a  useful  electrical  power  output  from 
the  external  circuit. 

The  above  scheme  accomplishes  the  primary  conversion  of  ther¬ 
modynamic  energy  directly  into  electrical  power  without  the  use 
of  any  major  rotating  or  reciprocating  components  such  as  large 
turbines  or  electrical  generators.  Nevertheless,  there  still 
exists  a  requirement  for  a  small  pump  to  recirculate  the  condensate 
and  a  fan  to  recirculate  the  carrier  fluid.  A  condenser  and  a 
boiler  are  also  normally  a  part  of  the  complete  system,  but  these 
various  auxiliaries  are  not  considered  here  in  any  detail  as  our 
present  concern  is  primarily  with  the  nozzle  in  which  the  basic 
thermo-electric  power  conversion  occurs. 

The  above  concepts  are  well  known  and  can  be  found  in  the 
technical  literature.  See  References  1,  2  and  3,  for  example. 

They  are  reviewed  in  this  introduction  merely  to  provide  a  proper 
background  for  our  result  which  is  quite  specific  and  which  is  a 
consequence  of  the  following  analysis. 

The  performance  that  can  be  achieved  by  an  EHD  device  is 
limited  by,  among  other  things,  the  maximum  electrical  field 
strength  that  can  be  sustained  at  the  most  critical  point  in  the 
field  without  inducing  electrical  breakdown,  through  a  spark 
discharge,  of  the  carrier  medium.  The  electrical  field  that  is 
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present  in  most  EHD  devices  has  its  highest  value  at  the  point 
of  injection.  See  References  1-4,  for  example.  In  this  regard 
there  is  experimental  evidence  that  breakdown  field  strength  in 
any  region  varies  linearly  with  the  quantity  Rp  where  R  is 
the  gas  constant  and  where  p  is  the  fluid  density  in  the  region 
considered.  This  experimentally  established  limit  depends  of 
course  on  the  type  of  gas/aerosol  combination  that  is  involved. 

But  the  presence  of  a  maximum  field  at  the  beginning  of  the 
EHD  conversion  channel  limits  performance  severely  because  the 
rest  of  the  channel  must  operate  below  maximum  capacity  or,  in 
other  words,  the  electric  pressure  is  at  its  allowable  limit  in 
only  a  small  fraction  of  the  conversion  channel. 

The  above  considerations  suggest  that  an  optimum  design 
would  be  one  in  which  the  local  field  strength  is  everywhere  uni¬ 
formly  close  to  the  critical  limit.  Such  a  design  would  achieve 
greater  electric  work  output  per  unit  mass  of  fluid  than  would 
any  other.  The  essence  of  this  paper  is  the  recognition  that  a 
uniform  maximum  loading  through  the  channel  is  the  optimum  load¬ 
ing  and  that  this  condition  can  indeed  be  achieved  at  least  to 
an  acceptable  approximation.  A  secondary  aim  of  this  paper  is  - 
to  derive  the  characteristic  geometrical  features  of  the  optimum 
design  and  certain  associated  performance  parameters  and  limits. 

The  geometrical,  electrical,  and  thermodynamic  features  of 
such  a  power  generator  are  governed  by  various  physical  laws  of 
which  one  of  the  most  significant  in  the  present  context  is 
Poisson's  equation  as  it  applies  to  an  electric  field.  In  its 
full  generality,  Poisson's  equation  is  three  dimensional  in 


nature;  by  restricting  attention  to  configurations  having  polar 
symmetry  we  can  simplify  this  to  a  two-dimensional  form.  More¬ 
over,  by  further  restricting  the  application  to  configurations 
whose  largest  radial  dimensions  are  small  in  comparisons  with 
their  axial  length,  we  may  finally  simplify  Poisson's  equation  to 
a  one-dimensional  approximation.  Despite  certain  limitations  in 
accuracy,  this  one-dimensional  version  can  be  more  enlightening 
than  the  more  elaborate  two-dimensional  analysis  because  it  shows 
basic  trends  so  much  more  simply  and  clearly.  It  is  this  one¬ 
dimensional  approximation  on  which  the  present  paper  is  based. 

The  analysis  must  deal  not  only  with  Poisson's  equation  for 
electrical  fields,  but  also  with  the  laws  of  fluid  flow.  It  is 
consistent  with  the  foregoing  commentary  to  write  these  laws  also 
in  their  common  one-dimensional  forms.  At  this  level  of  ideal¬ 
ization  it  is  also  appropriate  to  treat  the  flow  through  the 
channel  as  isen tropic.  We  also  consider  the  fluid  as  a  perfect 
gas  for  which  y  ,  the  ratio  of  specific  heats,  is  constant. 

One  other  special  assumption  is  involved.  This  is  based 
on  the  fact  that  the  drop  in  total  temperature  through  the 
channel  is  normally  very  small  in  comparison  with  the  absolute 
total  temperature  Tq  at  the  nozzle  inlet,  station  0.  Hence  it 
ia'  permissible  for  certain  analytical  purposes  to  neglect  this 
temperature  drop  and  to  approximate  the  flow  through  the  channel 
as  an  isentropic  flow  of  constant  total  temperature.  Calculations 
to  be  presented  later  show  that  the  actual  change  of  total  tem¬ 
perature  amounts  to  less  than  one  percent  for  a  typical  case. 
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Future  progress  in  EHD  power  generation  might  increase 
electric  power  output  by  an  order  of  magnitude  over  levels  that 
seem  currently  feasible.  In  that  case  the  drop  in  total  tempe¬ 
rature  through  an  EHD  nozzle  could  become  significant.  Fortunate¬ 
ly,  the  present  analysis  can  be  revised  and  generalized  to  account 
rigorously  for  such  variation  in  total  temperature  when  and  if 
necessary.  Under  present  circumstances ,  however,  it  is  not 
warranted  to  complicate  the  analysis  to  include  this  refinement 
since  at  low  power  outputs  it  has  no  appreciable  effect  on  the 
final  calculated  results. 

2.  Analysis 

It  is  convenient  for  the  purposes  of  the  present  analysis  to 
introduce  a  parameter  a  which  is  used  to  distinguish  between 
the  two  distinct  cases  of  positively  and  negatively  charged  par¬ 
ticles.  Specifically,  we  set  o  =  +  1  for  the  case  of  positively 

charged  particles,  and  we  set  a  -  -  1  for  the  case  of  negatively 

charged  particles.  It  is  then  appropriate  to  denote  the  electrical 
charge  per  unit  mass  by  means  of  the  product  cq  where  q  has 
units  of  coulombs/kg  and  is  always  positive  by  definition  while 
a  is  a  dimensionless  factor  as  defined  previously. 

Using  the  above  notation  we  may  write  the  following  two 
expressions  for  the  electric  current  i  and  the  mass  flow  rate 
ih  through  the  channel,  namely, 

i  »  aqpAv  »  constant  along  channel  (2.1) 

iti  =■  pAv  =  constant  along  channel  (2.2) 
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where 


i  =  electric  current,  amps 

ih  =  mass  flow  rate  through  channel,  kg/sec 

q  =  electric  charge  per  unit  mass,  positive  by 

definition,  coulomb/kg 

p  =  fluid  density,  kg/m3  (variable) 

A  =  cross-sectional  area  of  channel  at  any 

2 

given  streamwise  station,  m  (variable) 
v  =  mean  axial  velocity  at  any  given  streamwise 
station,  m/sec  (variable) 

Note  that  Eq.  (2.1)  is  based  on  the  previously  stated 
assumption  that  size  of  the  fluid  particles  and  the  charge  on 
each  particle  can  be  so  regulated  that  the  mobility  of  the  par¬ 
ticle  with  respect  to  the  surrounding  medium  is  essentially 
negligible . 

By  dividing  Eq.  (2.1)  by  Eq.  (2.2)  we  also  find  that 

=  aq  =  constant  along  channel  (2.3) 

The  one-dimensional  version  of  Poisson's  equation  can  now 
be  written  in  the  form 


(§H  - 


(2.4) 


where  z  =  axial  coordinate,  m 

0  =  electrical  potential,  volts 
£  =  eQ  =  electrical  permittivity  of  free  space 


=  8.854  x  10 


-2 


farad/m 
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The  local  electrical  field  strength  at  an  arbitray  axial 
station  z  is  defined  as  usual 

E  -  -  (i)  <2-s> 

where  E  =  field  strength,  volts/m 

As  will  be  seen  here,  an  optimum  EHD  generator  has  an  elec¬ 
tric  field  which  does  not  change  sign  from  inlet  to  outlet. 

Under  these  circumstances  the  electric  field  is  always  negative 
for  positive  space  charged  and  positive  for  negative  space  charge. 
Hence  we  may  write 

E  =  -  ct  |E  i  (2.6) 

It  has  been  established  by  experiment  that  over  a  broad  (but 
not  unlimited)  range,  the  field  strength  at  breakdown  is  well 
represented  by  the  simple  linear  law  (see  Reference  2) . 

iebI  -  co  +  S  R  »  <2-7’ 

where  C  and  C  are  characteristic  constants  of  the  medium. 

O  B 

R  is  the  gas  constant  in  units  of  Joule/kg°K  .  Experimental 
measurements  show  that  constants  Cq  and  CQ  happen  to  have  the 
same  numerical  values  for  both  air  and  steam.  However,  we  defer 
to  later  pages  of  this  paper  any  reference  to  actual  experimental 
values.  Of  course,  the  present  analysis  is  restricted  to  that 
specific  range  of  conditions  for  which  Eq.  (2.7)  is  in  fact  a 
valid  approximation. 

If  we  now  impose  the  constraint  that  the  local  field  strength 
is  everywhere  just  equal  to  its  critical  value  at  impending  break¬ 
down,  we  readily  find  from  the  last  three  equations  that 


(*) 


a(Co  +  CB  R  p) 


(2.8) 


Now  differentiating  this  gives 


/  d2<j>\ 
ydz2  f 


3  CB  R§ 


(2.9) 


Upon  equating  the  right  sides  of  Eqs.  (2.4)  and  (2.9)  and 
simplifying,  we  obtain  the  important  relation 


1  dp 
p"  dz 


■  CB  R 


(2.10) 


The  reciprocal  of  the  quantity  on  the  right  side  of  Eq.  (2.10) 
(momentarily  disregarding  the  negative  sign)  now  identifies  a 
significant  characteristic  length;  let  us  denote  it  by  symbol  X  . 


e  C  R 

X  =  -  (meters) 


(2.11) 


It  is  very  instructive  to  rewrite  Eq.  (2.11)  in  the  follow¬ 
ing  alternative  form 


Xq  =  X  . _  q  =  e  C_  R  =  constant 
n  min  ^max  B 


(2.12) 


Notice  that  the  quantity  e  C  R  is  a  characteristic  pro- 

O 

perty  of  the  medium.  Once  the  medium  is  chosen  the  designer  has 
no  further  control  of  the  value  of  this  constant.  Eq.  (2.12) 
now  tells  us  that  in  order  to  maintain  the  electrical  loading  at 
incipient  breakdown,  the  product  Xq  must  remain  constant.  Thus 
any  increase  in  q  must  be  accompanied  by  a  corresponding  decrease 
in  X  ,  or  vice  versa. 
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Experience  up  to  now  indicates,  however,  that  there  exists 

some  practical  upper  limit  <3max  to  the  value  of  q  that  can 

actually  be  achieved  for  any  specific  type  of  design,  and  hence 

there  exists  some  corresponding  lower  limit  X  .  on  the  accom- 

min 

panying  longitudinal  characteristic  length.  Of  course,  a  value 
of  q  less  than  qmax  can  always  be  employed  if  necessary  in 
which  case  the  corresponding  value  of  X  will  be  greater  than 

X  .  . 

min 

The  situation  is  shown  schematically  in  Fig.  2.1  in  the  form 
of  a  log-log  plot.  The  solid  line  represents  the  locus  of  points 
all  of  which  produce  incipient  breakdown.  Incidentally,  it  will 
also  be  shown  later  that  all  of  these  points  correspond  to  the 
same  theoretical  power  output  from  the  generator.  Note  that  the 
line  continues  on  indefinitely  toward  decreasing  values  of  q 
and  increasing  values  of  X  .  At  the  right  end,  however,  it  ter¬ 
minates  at  the  point  corresponding  to  q  and  X  ,  . 

"max  min 

Some  typical  values  of  these  various  characteristic  constants 
based  on  published  data  are  summarized  in  Table  2.1  for  reference 
purposes . 

Upon  substituting  X  into  Eq.  (2.10),  we  can  easily  inte¬ 
grate  the  simple  expression  that  is  thereby  obtained.  The  result 


X  (log  scale) 


Fig.  2.1  Design  Points  and  Design  Limits 

for  an  Optimum  Slender  EHD  Nozzle 


Table  2.1  Typical  EHD  Properties 


Source  of  Data 

'Ref.  1 

Ref.  2 

Ref.  3 

Medium 

Water 

droplets 

in  steam 

Water 

drop 16 ts 

in  air 

Water 

droplets 

in  air 

Quantity 

Units 

C 

o 

volts/m 

8.63  x  105 

8.63  x  105 

CB 

m2  °K/C 

9.49  x  103 

9.49  x  103 

9.49  x  103 

e  CB  R  =  Xq 

C  m/kg 

3.9  x  10"5 

2.4  x  10-5 

2.4  x  10"5 

^max 

C/kg 

2.5  x  10"3 

3  x  10~2 

3  x  10"2 

A  . 

min 

m 

1.5  x  10"2 

8  x  10'4 

8  x  10*4 

This  important  result  is  the  key  to  all  remaining  details  of 
the  analysis. 

Eq.  (2.13)  reveals  that  the  density  continues  to  drop  mono- 
tonically  with  increasing  distance  downstream.  Thus  the  flow  cor¬ 
responds  to  the  known  isentropic  flow  through  a  converging-diverging 
nozzle  for  which  the  walls  are  so  designed  that  the  density  drops 
exponentially  as  specified  by  Eq.  (2.13).  Let  subscript  1  denote 
the  throat  of  the  nozzle.  Then  the  following  relations  are  known 
to  apply 


1 


op 

I- 

1  +  Y"~  1  M2 

-< 

I 

*-* 

. 

(2.14) 
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-  1 
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+  T 


(Y  +  1) 

TTy  -■  IT 


l 

M 


(2.15) 


where  y  *  ratio  of  specific  heats 
M  »  Mach  number  at  station  z 
r  »  nozzle-  radius  at  station  z 
r^  «  nozzle  radius  at  throat 
It  is  readily  deduced  from  Eqs.  (2.13)  and  2.14)  that 


Cy  “=■  TT  £n 


(2.16) 


Eqs.  (2.15)  and  (2.16)  now  constitute  a  pair  of  parametric 

equations  for  (— )  as  a  function  of  (y)  with  Mach  number  M 
rl  A 

as  parameter.  These  equations  define  the  shape  of  the  optimum 
nozzle  having  uniform  loading  just  below  breakdown  all  along  its 

JT  Z 

length.  A  tabulation  and  a  curve  of  ( — )  versus  (y)  are 

rl  A 

given  in  the  next  section,  using  X  =  1.3  .  This  value  of  X  is 
representative  for  steam.  See  Table  4.1  and  Fig.  4.1. 

At  the  throat  of  the  nozzle  M  =  =  1  .  Then  Eq.  (2.16) 

gives 


For  y  *  1.3 


this  yields 


(2.17) 


0.466 


(2.18) 
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Let  z^  denote  the  axial  station  which  locates  the  injector 
that  introduces  unipolar  electrical  charges  into  the  flow.  We 
do  not  stipulate  in  advance  any  particular  fixed  location  for  the 
injector  because  the  nature  of  the  charged  particles  may  require 
one  location  or  another  depending  on  certain  supersaturation  con¬ 
ditions  that  might  be  needed  for  the  formation  of  the  aerosol 
droplets  (see  Reference  1) .  Moreover,  the  converging-diverging 
nozzle  does  not  have  to  be  of  the  shape  specified  in  this  analysis 
upstream  of  the  injector  location. 

Of  great  importance  is  the  theoretical  gross  power  output  of 
such  an  optimum  device.  To  determine  this  we  must  first  find  the 
current  from  Eq.  (2.1).  It  is  advantageous  to  evaluate  the  terms 
on  the  right  side  of  Eq.  (2.1)  at  the  nozzle  throat,  station  1. 
Thus, 

i  =  crqPjA^  (2.19) 

Now  letting  aQ  and  a^  denote  sonic  velocity  at  stations 
0  and  1,  respectively,  we  have 


Upon  combining  the  above  two  relations  we  finally  obtain  the 
current  in  the  form 


i 


aqpoaoAi 


(2.21) 
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Note,  however,  that  from  the  previous  analysis  of  the  isen- 
tropic  nozzle  flow,  it  follows  that 


(2.22) 


Moreover,  it  is  useful  to  eliminate  q  from  Eq.  (2.21)  by 
means  of  Eq.  (2.12)  and  to  rearrange  the  result  in  dimensionless 
form.  In  this  way  we  finally  obtain  the  dimensionless  current  in 
the  form 


/  Xi  ) 

•  " 

2 

V  k  B°aoAl )  " 

r  +  i 

.  a 

(y  + 


1) 

t r 


(2.23) 


The  next  step  is  to  find  the  change  in  electric  potential 
between  stations  i  and  2.  This  is  accomplished  by  integrating 
Eq.  (2.8)  between  these  limits.  We  again  rearrange  the  result 
in  dimensionless  form.  The  result  is 


(♦2  -  V  _ 

If  V  W»2  -  zt\  . 

■  "  A  X  -I 

A  _  a 

»o  S  R  X 

e  e 

Finally,  the  theoretical  power  output  follows  from  ttye  rela¬ 
tion 

=  i  ( <t>7  -  <fi.)  (2.25) 

e  *  i 

where  P@  *  theoretical  gross  electrical  power,  watts 
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The  desired  power  can  be  obtained  in  dimensionless  form  by 
direct  multiplication  of  Eqs.  (2.23)  and  (2.24).  After  some  minor 
rearrangement  we  obtain 


pe 

w 


-  e 


2 

Y  +  1 


(2.26) 


This  result  is  fundamental .  It  shows  just  how  the  theore¬ 
tical  gross  power  output  of  the  ideal  nozzle  depends  on  the  various 
geometrical  and  physical  parameters  of  the  problem. 

It  is  also  instructive  to  observe  from  Eqs.  (2.23),  (2.24) 
and  (2.26)  how  current,  voltage  and  power  vary  with  the  character¬ 
istic  length  parameter  X  of  the  nozzle.  In  other  words  we  hold 
nozzle  shape  constant  and  simply  vary  the  absolute  size.  Thus 
the  quantities  z^/X  and  z2/X  remain  constant.  Then  over  the 
range  X  >  X 


i  ~  X-1  (2.27) 

(<fr2  -  0i)  -  X  (2.28) 

Pe~  X°  =  independent  of  X  (2.29) 

Hence  current  is  inversely  proportional  to  size,  voltage 
is  directly  proportional  to  size,  and  power  is  independent  of 
size'.  This  comes  about,  of  course,  because  of  the  limitations 
imposed  by  the  electrical  breakdown  phenomenon*  see  also  Ref.  2. 
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Voltages  tend  to  be  inconveniently  high  in  EHD  devices. 

The  above  scaling  rules  suggest  that  this  problem  can  be  alle¬ 
viated  by  making  EHD  devices  small,  with  A  as  close  to  A  . „ 

nan 

as  feasible.  To  the  extent  that  it  can  be  practically  achieved, 
this  method  has  the  advantage  that  it  decreases  the  size  of  the 
device  without  any  corresponding  decrease  in  its  power  output/ 
Incidentally,  it  should  not  be  overlooked  that  Eq.  (2.26), 
being  based  on  a  purely  one-dimensional  approximation,  has  high 
accuracy  only  for  sufficiently  low  values  of  the  ratio  T2^z2  ' 
3.  Optimization 

It  is  convenient  to  simplify  Eq.  (2.26)  by  introducing  the 
following  auxiliary  notation.  Let 


(3.1) 


'  -  x 


and  let  Eq.  (2.26)  be  rewritten  in  the  form 


(3.2) 


_  ^  ^  £  _ 
^  R  po  ao  *1' 


•)  *  Fe =  [°^2  -  ci> +  e  ?1  - e  '2  (rnr) 


(y  +  D 
2(Y  -  1) 


(3.3) 


Another  quern tity  of  fundeunental  importance  is  the  exit 
kinetic  power  of  the  jet.  This  is  defined  as  follows 


Pj  ’  [»2  A2  V2] 


(3.4) 


By  using  the  various  well  known  relations  for  isentropic 
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compressible  flow  in  a  channel,  we  can  readily  reduce  this  to 
the  following  dimensionless  form 


(3. 


where  function  F^ (Mj)  is  as  follows: 
For  M2  <_  1 


M. 


(3y  -  1) 
27y  — ■) 


w 


^  V«2: 


(3. 


For  M2  >_  1 


F  .  (M,) 
3  2 


(Y  +  1) 

?rf  - ' rr 


- 1 


•l  Y  —  1  ..2 

1  +  — J—  M2 


(3. 


At  M2  =r  1  -these  two  expressions  both  give  the  result 

{ 3y  -  1) 

2(y-"-  i) 

P3I1>-i[y4t]  <3 

It  will  also  be  useful  later  to  have  the  value  of  F ^ 
in  the  limit  as  M2  increases  without  bound.  From  Eq.  (3.7) 
we  may  deduce  that  this  is 


(Y  +  1) 


Consider  an  EHD  power  generator  which  operates  as  a  closed 


system.  The  fluid  which  discharges  from  the  nozzle  exit  must 
be  recirculated  back  to  the  nozzle  inlet.  This  involves  dece¬ 
lerating  the  flow,  cooling  the  fluid,  condensing  and  separating 
the  liquid  droplets,  and  so  on.  It  involves  restoring  the 
pressure  of  the  fluid  by  means  of  a  suitable  fan,  compressor  and/ 
or  pump,  reheating  the  liquid  and  completing  the  thermodynamic 
cycle.  These  processes  require  power.  The  basic  compressor  and/ 
or  pump  power  requirement  is  inherent  in  the  thermodynamic  cycle 
itself  and  cannot  be  eliminated.  Additional  power  demands  arise 
from  the  various  irreversibilities  that  inevitably  occur  in  any 
real  physical  system;  this  additional  power  cannot  be  completely 
eliminated  but  it  can  be  minimized  by  careful  engineering  design. 
Fortunately,  some  of  the  power  required  can  actually  be  recovered 
from  the  exit  kinetic  power  of  the  jet  itself.  Inevitably,  how¬ 
ever,  there  will  be  a  net  power  demand,  call  it  5  P^  ,  which  must 
be  subtracted  from  the  gross  electric  power  output  Pg  .  Conse¬ 
quently,  the  net  useful  electrical  power  becomes 


P 


e  net 


P 

e 


5  P. 
1 


(3.10) 


Yet  another  jet  power  term  is  very  useful.  This  is  the  jet 
power  developed  at  the  throat  of  the  nozzle,  call  it  Pj (1)  . 
This  makes  a  useful  dimensional  reference  parameter  which  may  be 
used  to  normalize  the  net  power  output  Pg  net  . 

This  reasoning  leads  us  to  define  the  following  useful  net 
power  coefficient,  namely. 
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X 


(3.11) 


Pe  net 

T-TTT 


-  5 

T~m 


Upon  substituting  into  this  definition  the  expressions  for 
Pg  ,  Pj  and  Pj(l)  as  previously  developed,  and  upon  simpli¬ 
fying,  we  can  reduce  the  result  to  the  form 


X 


1 

fTTTT 


8  F 


(3.12) 


(3.13) 


It  is  therefore  clear  that  we  should  maximize  inlet  stagna¬ 
tion  pressure  pQ  while  keeping  Tq  no  higher  than  necessary  to 
obtain  proper  droplet  size  in  the  medium.  Moreover  it  might  be 
possible  to  eliminate  this  restriction  on  Tq  by  injecting  solid 
or  liquid  charged  particles  rather  than  condensing  droplets.  In 
this  case  Tq  could  be  sharply  reduced  thus  greatly  incresing 
the  value  of  the  crucial  parameter  8  . 

Eq.  (3.12)  for  the  power  coefficient  x  is  the  key  result 
of  this  section.  The  optimization  of  an  EHD  generator  can  be 
based  in  part  on  an  effort  to  maximize  this  quantity.  In  this 
connection,  however,  another  quantity  of  interest  is  a  power 
ratio  n  which  we  define  as  follows : 


(3.14) 
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Notice  that  the  quantity  n  as  defined  above  does  not 
allude  to  the  net  heat  input  to  the  cycle.  Hence  it  should  not 
be  confused  with  overall  thermal  efficiency  which  will  always 
be  some  very  much  lower  figure. 

The  various  power  terms  that  occur  in  Eq.  (3.12)  can  also 
be  used  to  determine  the  ratio  of  ATq  ,  the  drop  in  total  tem¬ 
perature  through  the  nozzle,  Tq  ,  the  inlet  total  temperature. 
This  works  out  to  be  simply 


AT 


B  Fe 

F  j  (  ®)" 


(3.15) 


Recall  that  the  present  analysis  is  restricted  to  cases  in 
which  AT  /Tq  is  small  compared  with  unity.  Eq.  (3.15)  may  be 
used  to  check  whether  this  restriction  is  in  fact  satisfied  in 
axny  particular  case. 

It  is  evident  from  Eq.  (3.12)  that  the  value  of  the  net 
power  coefficient  x  that  can  be  attained  is  governed  by  the 
two  key  parameters  8  and  <5  .  The  coefficient  6  controls 
the  maximum  gross  power  that  can  be  generated  without  electrical 
breakdown,  and  the  factors  that  govern  8  in  turn  are  shown  in 
Eq.  (3.13).  Coefficient  6  describes  the  power  loss  inherent 
in  the  thermodynamic  cycle  itself,  as  augmented  by  the  various 
irreversibilities  that  must  inevitably  occur  in  such  a  system. 

4.  Summary  of  Key  Relations  and  Results 

The  maximum  specific  charge  q  that  can  be  employed  in  an 
ideal,  slender  EHD  nozzle  without  causing  electrical  breakdown 
is  governed  by  the  relation 
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=  Amin  \ax  -  £  CB  R  *  constant 


The  proper  geometric  form  for  an  optimum  EHD  no2zle  is 
defined  by  the  following  pair  of  parametric  equations 

(Y  +  1) 

,  V  r,  .  A  -  1  ,.2t  “T  *  15 


/  v  r-i  .  A  —  1  ..2~i 

/r\  1  +  ~T~  M 

\ri)~  I  rip 
c  =  |l j=  TT^TT  *n 


in  i  +  r  :  .l  m2 


The  value  of  $2  corresponding  to  any  chosen  value  of  exit 
Mach  number  may  be  calculated  from  Eq.  (4.3) .  Then  we  set 


a  . 


B  oi 


and  find 


*u2  -  cA 

. 


-  Ci  -  c. 


)  +  e  1  - 


(Y  +  1) 

Ty  -  ~ir 


where  denotes  the  location  of  the  injector. 

The  following  functions  and  constants  may  next  be  found. 
For  M2  <_  1 

(3y  -  1) 


M  3  T  /  V  I  2(Y 

’^H^R 


From  >  1 


Finally 


(4.7) 

(4.8) 

(4.9) 

(4.10) 

(4.11) 

(4.12) 


AT  8  F 

T  =  p  (co)  (4.13) 

o  j 

Values  of  x  /  n  and  ATq/To  versus  M2  are  listed  in 
Table  4.2  for  a  typical  case  which  is  defined  by  the  specified 
values  of  y  ,  ^  ,  a  ,  8  and  5  .  The  value  8  =  0.1  corres¬ 
ponds  to  air  or  steam  at  about  250  atmospheres;  this  example 
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illustrates  the  extreme  measures  that  are  required  to  lift  per¬ 
formance  to  levels  which  are  even  marginally  acceptable.  The 
assumed  value  5  *  0.1  implies  thac  pump  work  and  all  fluid 
losses  amount  to  only  10  percent  of  the  kinetic  power  of  the  jet; 
this  is  a  very  optimistic  hypothetical  example,  but  it  suffices 
to  show  certain  trends.  Reference  (5)  gives  more  detailed  data 
on  one  particular  type  of  loss  that  may  be  involved,  namely, 
losses  in  pipe  or  channel  bends. 

Table  4.2  shows  clearly  the  general  effects  of  varying  exit 
Mach  number  M2  .  It  also  shows  the  trade-off  between  perfor¬ 
mance  parameters  x  and  n  ;  generally  speaking  an  increase  in 
either  of  these  can  only  be  obtained  at  the  expense  of  a  decrease 
in  the  other.  Note  that  excessively  high  values  of  M2  lead  to 
negative  values  of  outputs 

The  last  column  of  Table  4.2  reveals  that,  over  the  range 

of  greatest  interest,  the  value  of  AT  /T  is  smaller  than  0.01. 

o  o 

This  confirms  that  one  of  the  basic  restrictions  of  the  theoretical 
analysis  is  indeed  satisfied  in  this  case. 

The  optimum  EHD  device  disclosed  herein  is  useful  not  only 
as  a  possible  real  physical  system  but  also  as  a  hypothetical 
standard  against  which  the  performance  of  any  other  real  or  hypo¬ 
thetical  EHD  system  may  be  compared. 

The  analysis  of  this  section  could  be  greatly  improved  by 
treating  the  optimum  EHD  channel  as  just  one  element  of  a  com¬ 
plete  thermodynamic  cycle.  By  working  with  such  a  complete  cycle 
in  detail  the  true  value  of  parameter  6  ,  instead  of  being 
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TABLE  4.1 


OPTIMUM 

SLENDER 

EHD  NOZZLE 

M 

Y  *  1.3 
z/X 

).  2° 

0.  010  3 

1. 7303 

0.4  0 

0.  0  701 
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1.6  0 

1.  03  33 
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1.00 

1. 3203 

1. 2132 

2.00 

1.5557 

1. 3315 

2.20 

1.3104 
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2.4  3 

2.0757 

1 .  -5  2  J  0 

2.00 

2.3337 

1.3153 

2.00 

2 . 5  0 16 

2. 0230 

3.00 

2. 3431 

2.2715 

3.20 

3.1020 

2  .  5  4  51 

3.4  0 

3  .  3  32-3 

2.3513 

3 . 5  0 

3.5)02 

3.1)33 

3.30 

3 .  3  415 

3.5734 

4.  00 

4. 0733 

3.0)30 

4.20 

4.3121 

4.4551 

4.40 

4. 5400 

4.0521 

4 . 3  0 

4.7523 

5.5155 

4. 3  0 

4.3203 

3. 1214 

5.00 

5.1333 

5.7701 
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, «  ».r'  ^ 
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TABLE  4.2  PERFORMANCE  OF  OPTIMUM  EHD  GENERATOR 


Y  =  1.3  =  0  a  =  0  8  =  0.1  6  =  0.1 


M2 

C2 

X 

n 

‘VTn 

0.05 

0.0012 

0.0003 

0.9146 

0.0000 

0.10 

0.0050 

0.0010 

0.3297 

0,0001 

0.15 

0.0112 

0.0019 

0 . 7459 

0 . 0003 

0.20 

0,0199 

0.0030 

0 , 6636 

0*0 0 0  6 

0.25 

0.0311 

0,0041 

0.5333 

0.0009 

0.30 

0.0447 

0.0051 

0 , 5054 

0 . 00 1 3 

0.35 

0.0607 

0 . 0058 

0.4304 

0.0013 

0.40 

0.0791 

0 . 0063 

0.3585 

0.0023 

0.45 

0.0997 

0 . 0063 

0.2903 

0.0023 

0.50 

0.1227 

0.0060 

0 . 2260 

0.0035 

0 . 55 

0.1479 

0,0052 

0 . 1653 

0.0041 

0 . 60 

0. 1753 

0.0041 

0.1099 

0.0043 

0 . 65 

0.2043 

0 . 0025 

0,0587 

0  » 0056 

0.70 

0.2364 

0 , 0006 

0,0121 

0.0063 

0 . 75 

0.2700 

-0.0016 

-0.0297 

0.0071 

0.80 

0 . 3056 

-0.0040 

-0 » 0668 

0 . 0079 

0.35 

0 . 3430 

-0.0066 

-0.0990 

0 , 0087 

0.90 

0 . 3822 

-0.0092 

-0.1265 

0.0095 

0.95 

0.4232 

-0.0113 

-0.1493 

0.0104 

1.00 

0.4659 

-0.0143 

-0.1675 

0,0112 

>> 
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merely  assumed,  could  be  actually  calculate  for  each  case.  Also 
instead  of  working  with  relative  indices  of  performance  like  x 
and  n  ,  we  could  calculate  actual  overall  thermodynamic  effi¬ 
ciency.  Moreover,  the  effect  on  overall  thermodynamic  efficiency 
of  varying  various  characteristics  of  the  thermodynamic  cycle 
could  then  be  systematically  evaluated.  Time  limitations  have  pre¬ 
cluded  such  a  more  ambitious  and  informative  effort.  The  present 
report  concentrates  instead  primarily  on  the  concept  of  the  opti¬ 
mum  EHD  channel  itself. 

5.  Comparison  of  Optimum  Channel  with  Uniform  Channel 

In  EHD  power  generation,  a  uniform  channel  has  the  obvious 
merit  of  simplicity.  In  this  case  the  gas  density  remains  prac¬ 
tically  constant.  However,  the  electrical  field  is  always  a  maxi¬ 
mum  at  the  entrance  to  the  working  section,  and  gradually  decreases 
thereafter.  Maximum  power  output  is  obtained  if  the  length  be  so 
chosen  that  the  electrical  field  is  zero  at  the  exit  from  the 
working  section. 

It  is  well  known  that  the  electrical  power  output  of  such  a 
uniform  channel  is 

Pe  -  y  e  EQ2  AV  (watts)  (5.1) 

where 

Eg  »  breakdown  strength  of  medium,  volts/m 

2 

A  »  area  of  uniform  channel,  m 

V  *  constant  mean  velocity  in  channel,  m/sec 


For  the  purpose  of  the  present  section  we  divide  Pg  by 
the  mass  flow  rate  in  thus  fixing  the  electrical  work  per  unit 
mass  of  fluid  (Joule/kg) .  It  is  also  convenient  to  nondimen- 
sionalize  the  resulting  quantity  by  dividing  through  by  a*  /2 
(Joule/kg)  where  a*  is  the  sonic  velocity  which  corresponds  to 
a  Mach  number  of  unity.  We  also  utilize  Eq.  (2.7)  and  make  use 
of  the  fact  that  the  constant  CQ  is  negligible  in  comparison 
with  the  quantity  CQRp  •  In  addition  we  employ  various  standard 
perfect  gas  relations.  In  this  way  we  reduce  Eq.  (5.1)  to  ob¬ 
tain  a  dimensionless  power  coefficient  C  as  follows: 

6 


This  shows  clearly  that  reducing  the  Mach  number  M  increases 
the  electric  work  per  unit  mass  of  fluid.  The  maximum  value  that 
can  be  achieved  in  this  way  is 


Now  let  us  consider  the  optimally  loaded  configuration. 

Eq.  (2.26)  provides  a  convenient  starting  point.  As  in  the  pre¬ 
vious  case  we  neglect  the  quantity  (C  /C_  R  p  )  thus  eliminating 
the  corresponding  term  from  the  equation.  The  maximum  power  is 
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The  result  can 


obtained  by  setting  z^/\  *  0  and  z ^ /*  *  ®  • 
then  be  nondimensionalized  in  much  the  same  way  as  for  the  pre¬ 
vious  case.  Notice  that  in  this  case  sonic  velocity  occurs  at 
the  throat  of  the  nozzle  so  that  a*  *  a1  .  The  required  result 
can  now  be  developed  as  follows 


The  bracketed  exponent  turns  out  to  be  zero  whereupon  the 
final  result  becomes  simply 


Inspection  of  Eqs.  (5.3)  and  (5.5)  now  reveals  that  the  op¬ 
timally  loaded  channel  produces  twice  as  much  electric  work  per 
unit  mass  of  fluid  as  does  the  uniform  channel.  These  formulas  - 
are  of  fundamental  importance  because  they  show  the  basic  per- 
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formance  limits  of  these  two  types  of  EHD  channels.  Substitu¬ 
tion  of  typical  numerical  values  into  these  formulas  will  also 
show  how  stringent  these  performance  limits  really  are. 

An  improvement  in  performance  by  a  factor  of  2  is  unfortu¬ 
nately  not  sufficient  in  itself  to  bring  EHD  power  generation  to 
practical  reality.  A  further  improvement  by  a  factor  of  somewhere 
between  5  and  50  is  probably  still  needed  to  achieve  that  goal. 
Moreover  a  factor  of  2  improvement  may  not  be  sufficiently  profit¬ 
able  under  the  circumstances  to  warrant  the  extra  geometrical  and 
other  complexities  involved.  Nevertheless,  the  two  performance 
limits  developed  above  have  considerable  theoretical  interest  and 
value. 

Incidentally,  recall  that  another  simple  but  important  result 
of  the  present  paper  is  its  identification  of  the  key  dimension¬ 
less  scaling  parameter  X  .  Eqs.  (2.27)  through  (2.29)  show  the 
basic  significance  of  this  quantity. 
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